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STRAIN-RATE  EQUATIONS  FOR  CALCULATION  OF  SECONDARY  CREEP 
DEFORIVIATION  OF  THICK-WALLED  TUBES  WITH 
INTERNAL  OR  EXTERNAL  PRESSURE 
by  Richard  E.  Morris 
Lewis  Research  Center 

SUMMARY 

The  design  of  tubular  components  of  nuclear  propulsion  systems  for  long  service 
life  at  high  temperature  and  pressure  requires  that  the  secondary  creep  deformation  be 
controlled  or  limited.  Examples  of  such  structures  include  heat  exchangers  and 
nuclear-reactor  fuel-pin  cladding. 

Heat  exchangers  must  operate  with  internal  pressure  at  high  temperature  for  long 
lifetimes.  An  important  requirement  is  the  limitation  of  the  creep  expansion  of  the 
thick-walled  heat-exchanger  tubing. 

In  the  case  of  nuclear-reactor  fuel-pin  cladding,  both  external  coolant  pressure 
and  the  development  of  internal  fission  gas  pressure  cause  creep  deformations  that 
must  be  limited. 

Strain-rate  equations  are  presented  for  use  in  these  design  problems.  The  stress 
and  strain-rate  relations  apply  to  thick-walled  tubes  for  both  internal  and  external 
pressure.  Both  temperature  and  pressure  were  assumed  to  be  constant.  Stress  and 
strain  rate  were  assumed  to  be  related  by  a  power  law  which  includes  a  function  of 
temperature  as  a  parameter.  Thus,  the  equations  are  applicable  to  the  stress  and 
strain-rate  behavior  of  tubes  over  a  range  of  temperatures. 


INTRODUCTION 

Tubular  components  of  nuclear  propulsion  systems  may  be  designed  to  operate 
under  extreme  conditions  of  temperature  and  pressure  such  that  creep  deformation 
occurs  continuously  during  operation. 


An  important  limitation  in  designing  heat  exchangers  for  high-temperature  opera¬ 
tion  is  the  creep  strength  of  heat-exchanger  materials.  Similarly,  the  creep  strength 
of  nuclear -reactor  fuel-pin  cladding  materials  poses  a  limitation  on  the  operating  tem¬ 
perature  of  the  fuel  pins  in  a  reactor. 

Helium-to-air  heat  exchangers  may  operate  with  an  internal  helium  pressure  of 
7  to  14  meganewtons  per  square  meter  at  temperatures  from  1000  to  1150  K  for  life¬ 
times  of  50  000  hours.  Under  these  conditions  the  heat-exchanger  tubes  undergo  con¬ 
tinuous  creep  deformation.  Rupture  of  heat-exchanger  tubes  must  be  avoided.  The 
amount  of  deformation  accumulated  during  the  operating  lifetime  of  the  tubing  must  be 
limited  by  design. 

Tubular  nuclear -reactor  fuel-pin  cladding  contains  the  nuclear  fuel.  Fission  gases 
are  generated  as  the  fuel  is  consumed.  Gases  collect  in  the  high-temperature  core 
void  of  a  fuel  pin.  The  gas  is  constrained  by  the  volume  of  void  provided  in  the  fuel- 
pin  design  so  that  pressure  develops  as  a  function  of  fuel  burnup.  Pressures  as  high  as 
25  meganewtons  per  square  meter  may  be  generated  during  the  life  of  a  pin. 

For  high-performance  reactors  the  peak  temperature  of  fuel-pin  cladding  may  be 
from  1100  to  1500  K.  In  this  temperature  range,  the  regularly  increasing  fission  gas 
pressure  will  cause  the  cladding  to  creep  during  the  10  000-hour  design  life  that  is 
typical  of  these  fuel  pins.  Radioactive  fission  gases  must  be  contained  within  the  fuel 
pin.  Rupture  must  be  avoided.  Total  creep  deformation  during  the  lifetime  of  the  fuel- 
pin  cladding  must  be  limited  by  design. 

The  coolant  in  the  gas-cooled  reactor  must  operate  under  pressure  sufficient  to 
obtain  the  heat-transfer  characteristics  needed  to  cool  the  fuel  pins.  For  heliiim,  a 
coolant  pressure  of  about  10  meganewtons  per  square  meter  is  needed.  This  pressure 
is  sufficient  to  cause  continuous  compressive  creep  of  the  fuel-pin  cladding  during  the 
first  part  of  the  life  of  the  reactor.  Compressive  creep  must  be  limited  to  avoid  col¬ 
lapse  of  fuel-pin  cladding.  The  fission  gas  pressure  gradually  increases  to  oppose  the 
coolant  pressure.  Fission  gas  pressure  is  greater  than  the  coolant  pressure  at  the  end 
of  the  life  of  the  reactor. 

Each  of  these  design  problems  involves  the  relation  between  stress  and  strain  rate 
in  thick- walled  tubes  undergoing  creep  from  either  external  or  internal  pressure. 

This  report  presents  stress  equations  and  strain-rate  equations  for  applications  to 
these  problems.  The  analysis  follows  the  method  presented  by  Bailey  (ref.  1)  which  is 
a  three-dimensional  stress  analysis  of  internally  pressurized  tubes  undergoing  second¬ 
ary  creep.  The  analysis  presented  in  this  report  changes  one  of  his  fundamental 
assumptions  so  that  the  strain-rate  equations  apply  to  a  range  of  operating  tempera¬ 
tures.  In  addition,  by  using  the  boundary  conditions  for  external  pressure  the  analysis 
yields  stress  equations  and  strain-rate  equations  for  tubes  under  external  pressure. 

All  of  the  equations  apply  only  to  the  case  of  constant  pressure  and  temperature. 
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SYMBOLS 


1  9  —n 

A  material  constant,  hr”‘‘'(N/m^)~'^ 

a  inside  radius,  cm 

B  material  constant,  hr”'*'  (N/m^)'** 

b  outside  radius,  cm 

C  constant  of  integration 

D  constant  of  integration 

AH  apparent  activation  energy,  J/(K)(mole) 

n  stress  exponent 

2 

p  pressure,  N/m 

R  gas  constant,  8.  3143  J/(K)(mole) 

r  radius,  cm 

T  absolute  temperature,  K 

u  radial  displacement,  cm 

•  _  1 

e  creep  rate,  hr 

e  equivalent  creep  rate,  hr"^ 

P  b/a 

2 

a  stress,  N/m 

—  2 
a  equivalent  stress,  N/m 

X  b/r 

Subscripts: 

i  internal 

0  external 

r  radial 

z  axial 

9  circumferential 


ANALYSIS 


This  analysis  follows  the  work  of  Bailey  for  the  three-dimensional  stress  analysis 
of  internally  pressurized  tubes  (ref.  1)  except  for  the  power  law  assxunption  that  re¬ 
lates  stress  and  strain  rate.  Bailey  assumed  a  power  law  e  =  Ao^.  In  this  form  both 
A  and  n  are  temperature  dependent.  Consequently,  strain-rate  relations  for  creep 
in  thick-walled  tubes  under  internal  pressure  based  on  this  power  law  are  applicable 
only  for  a  relatively  narrow  range  of  temperature  for  which  A  and  n  are  constant. 

The  power  law  assumed  in  this  analysis  contains  a  temperature  correction  term 
g-AH/RT  ygg  term  reduces  the  temperature  dependence  of  the  constants  in 

the  stress,  creep-strain-rate  relation.  With  the  temperature  correction  term  included, 
the  creep-strain-rate  relation  used  in  this  analysis  is 


This  power  law  has  been  used  by  several  investigators,  for  example,  Maag  and  Mattson 
(ref.  2).  Garofalo  (ref.  3)  indicates  that  this  strain-rate  relation  applies  to  secondary 
creep  at  constant  temperature  for  many  metals  and  alloys  at  low  stresses.  Long  life¬ 
times  with  low  creep  rates  at  high  temperature  Limit  the  use  of  many  metals  to  low 
stresses. 

The  use  of  the  power  law  with  the  temperature  correction  term  in  this  analysis  and 
adding  the  boundary  conditions  for  tubes  with  external  pressure  provide  strain-rate 
relations  for  thick-walled  tubes  loaded  with  static  internal  or  external  pressure  under¬ 
going  isothermal  secondary  creep  as  a  function  of  temperatures. 

The  assumptions  used  in  the  analysis  are  as  follows: 

(1)  The  tube  material  is  isotropic. 

(2)  Creep  associated  with  the  deformation  of  the  tubes  is  secondary  creep. 

(3)  Secondary  creep  strain  follows  the  von  Mises  flow  rule. 

(4)  Stress  and  strain  rate  are  assumed  to  be  related  by  the  empirical  power  law 

e”  =  (The  sign  is  selected  for  tensile  (+)  or  compressive  (-)  strain  rate. ) 

(5)  The  axial  creep  rate  is  assumed  to  be  zero.  Planes  perpendicular  to  the  axis 
of  the  tube  remain  plane.  The  creep  deformation  is  one  of  plane  strain. 

(6)  The  principal  axes  of  stress  and  creep  strain  rate  are  coincident  and  remain 
so  during  strain  up  to  moderate  amounts. 

(7)  The  principal  shear  strain  rates  are  proportional  to  the  principal  shear  stresses. 

(8)  Poisson's  ratio  is  assumed  to  be  1/2.  This  assumption  requires  that  the  sum 

of  principal  orthogonal  creep  rates  equal  zero,  that  is  ~ 
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INTERNAL  PRESSURE 


The  use  of  the  power  law  assumption  with  the  temperature  correction  factor  results 
in  strain-rate  equations  for  and  Correlation  of  the  equations  with  experi¬ 

mental  data  will  provide  values  for  the  material  constants.  The  resulting  equations 
are  more  general  than  Bailey's  strain-rate  equations,  in  that  the  equations  presented 
here  include  temperature  as  a  parameter. 

The  stress  equations  derived  by  Bailey  (ref.  1)  are  not  affected  by  this  change  in 
the  stress-strain  rate  power  law  assumption.  The  equations  for  the  stresses  in  a  tube 
undergoing  secondary  creep  at  constant  temperature  are  given  here  in  the  terminology 
of  this  report  for  completeness  (the  derivation  is  included  in  the  appendix); 


2/n 


p2/n.l  ^ 


1  _  2/n 


^0  = 


n 


,2/n 


CT. 


n 


2/n 


(1) 


(2) 


(3) 


(4) 


Equation  (4)  is  obtained  by  substituting  the  stress  equations  into  the  relation  for  the 
maximum  distortion  energy  stress  (equivalent  stress) 


a=-J-[( 

V2 


(a^.  -  ag)2  +  (ag  -  +  (a^ 


-  o^)\ 


1/2 


Use  of  the  modified  power  law  assumption  yields  the  following  strain-rate  equations: 


■■  2 


(5) 
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e  -  (6) 

^  2 

=  0  (assumption  (5))  (7) 

z 

?  =  (assumption  (4),  tensile  creep  rate)  (8) 

These  eight  equations  define  the  stresses  and  strain  rates  in  an  internally  pressur¬ 
ized  tube  at  constant  temperature  for  secondary  creep  conditions. 


External  Pressure 

The  derivation  for  stresses  in  tubes  from  internal  pressure  uses  the  boundary 
conditions 


=  -p  at  r  =  a 

Oj,  =  0  at  r  =  b 

Changing  the  boundary  conditions  to 

-  0  at  r  =  a 

=  -p  at  r  =  b 

we  have  the  condition  for  tubes  loaded  by  external  pressure.  Then  if  the  boundary  con¬ 
ditions  are  used  to  evaluate  the  constants  of  integration  in  the  derivation,  the  equations 
for  the  stresses  become 


2/n  _  2/n 

a  -  -  I- - ^ - p^ 

O  /»,  0 

-  1 


p2/n  - 

-  1 


(9) 


(10) 
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(11) 


p2/n-MU2/n 


n 


p2/n  -  1 


p2A.l  ° 


(12) 


Assiunption  (4)  is  modified  for  negative  strain  rates  assuming  that  compressive 
behavior  is  similar  to  tensile  behavior.  The  resulting  strain  rates  are 


e  (13) 

^  2 

2 


e_  =  0  (assumption  (5))  (15) 

?  =  (assumption  (4),  compressive  creep  rate)  (16) 

Equations  (13)  to  (15)  relate  the  stresses  and  strain  rates  to  functions  of  pressure, 
geometry,  temperature,  and  materials  properties  for  a  given  radius  in  the  wall  of  a 
tube.  The  sign  of  the  empirical  exponential  stress  function  is  assumed  to  be  posi¬ 
tive. 


DISCUSSION 

Equations  are  presented  for  the  stresses  and  strain  rates  in  isothermal  thick- 
walled  tubes  with  static  internal  or  external  pressure  that  are  applicable  over  a  range 
of  temperatures.  Tube  structures  that  are  designed  for  long  lifetimes  at  high  tempera¬ 
tures  may  undergo  creep  deformation  continuously.  Design  stresses  must  be  selected 
so  that  creep  rates  are  limited  to  very  low  secondary  creep  rates.  If  this  is  done,  ter¬ 
tiary  creep  deformation  will  be  avoided.  The  equations  presented  apply  only  for  tubes 
undergoing  secondary  creep  deformation. 

The  power  law  between  stress  and  strain  rate  used  in  developing  the  equations  con¬ 
tains  a  temperature  correction  term.  Thus,  the  equations  obtained  are  more  general 
than  previous  work  because  they  can  be  used  to  describe  strain-rate  data  for  a  range  of 
temperatures. 
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The  equations  presented  constitute  a  closed-form  solution  to  the  problem  of  iso¬ 
thermal  secondary  creep  in  thick-walled  tubes  loaded  with  either  internal  or  external 
pressure.  When  the  material  properties,  constants,  and  tube  geometry  are  substituted 
into  the  equations,  the  creep  strain  rates  and  stresses  are  completely  described 
throughout  the  wall  of  the  tube. 

Fuel-pin  cladding  and  heat-exchanger  tubing  are  generally  designed  for  thousands 
of  hours  of  operation.  The  principal  source  of  long-term  deformation  of  these  tubular 
structures  will  be  secondary  creep  from  pressures  acting  on  the  structures. 

Some  observations  can  be  made  by  comparing  the  equations  for  stress  and  strain 
rate  for  external  and  internal  pressure  loading.  The  equations  for  strain  rate  are 
similar  but  of  opposite  sign.  The  equations  for  stresses  ct^.,  a^,  and  are  different, 
while  the  equations  for  the  equivalent  stresses  are  the  same  but  of  opposite  sign.  If  the 
internal  pressure  and  the  external  pressure  acting  on  a  tube  are  equal,  the  net  equiva¬ 
lent  stress  will  be  zero  everywhere  in  the  tube  wall.  Consequently,  no  creep  will  occur. 

CONCLUDING  REMARKS 

Temperature-dependent  strain-rate  equations  are  presented  for  thick-walled  tubes 
undergoing  isothermal  secondary  creep  from  internal  pressure.  Temperature- 
dependent  stress  and  strain-rate  equations  are  provided  for  thick- walled  tubes  under¬ 
going  isothermal  secondary  creep  from  external  pressure. 

These  equations  are  applicable  to  the  analysis  of  stresses  and  long-term  isothermal 
creep  deformation  of  thick-walled  tubular  structures  loaded  by  static  internal  or  ex¬ 
ternal  pressures. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  June  10,  1971, 

126-15. 
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APPENDIX  -  DERIVATION  OF  STRESS  AND  STRAIN-RATE  EQUATIONS  FOR 
THICK-WALLED  TUBES  AT  CONSTANT  TEMPERATURE 
UNDERGOING  SECONDARY  CREEP 
FROM  INTERNAL  PRESSURE 


The  derivations  in  this  appendix  are  from  the  analytical  precedure  presented  by 
Bailey  in  reference  1.  The  assumption  that  e  =  0  is  in  agreement  with  experimental 
evidence.  The  assumption  implies  that  planes  perpendicular  to  the  axis  of  a  tube  re¬ 
main  plane.  This  is  the  condition  of  plane  strain. 

Force  equilibrium  for  cylindrical  geometry  in  the  absence  of  body  force  yields 


do-j. 

dr 


"  -  '^r 


The  strain  displacement  relations  are 


"0  = 


u 


"r  = 


dr 


Differentiating  with  respect  to  time  yields 


S  = 


dr 


(Al) 


(A2) 


(where  dots  indicate  d/dt).  The  compatibility  equations  for  this  case  reduct  to 


r 


dr 


(A3) 
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Creep  is  a  constant  volxime  process,  so  that 


Applying  the  assumption  that  e  =  0,  we  have 


Substituting  equation  (A4)  into  (A3)  yields 


e  2dr 


Upon  integration,  we  have 


€/,  =  Cr 


where  the  constant  C  is  to  be  evaluated. 

The  assumption  that  the  principal  shear  strain  rates  are  proportional  to  the  princi 
pal  shear  stresses  can  be  expressed  as  follows; 


~  ^r  _  ^r  ~  ^z  _  ^z  ~ 
^9-^r  '^z-^f 


Substituting  equation  (A4)  into  (A6)  yields 


An  expression  for  the  equivalent  stress  cr  is  given  by  the  distortion  energy  theory 
(also  known  as  the  von  Mises  yield  criterion): 


^  [(<^0  -  +  (""r  -  ^  ■  ""0^^] 
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Substituting  for  a_  from  equation  (A7)  gives 


cr  =  ^  (cTq  -  cr^) 


(A9) 


The  equivalent  strain  rate  e  corresponding  to  the  equivalent  stress  is  given  by 

1 1/2 


(AlO) 


Substituting  for  e^  =  0  and  =  -e^  reduces  equation  (AlO)  to 


^  2 


V3 


e 


(All) 


Make  the  assumption  that  experimental  strain-rate  data  can  be  represented  by 

?  =  (A12) 


Then 


B 


2 


ng-AH/RT 


(A12) 


and  from  equations  (A5)  and  (Al) 


(A13) 


Let 


D  = 


B 


|l/n 


11 


daj.  =  D 


dr 

^(n+2)/n 


Integration  yields 


1 

2/n 
r  ' 


+  F 


The  constants  of  integration  C  and  F  are  evaluated  using  the  boundary  conditions . 
At  r  =  b,  =0 


+  F  =  0 


Dn 

2b2A 


a. 


Dn/  1  1  \ 

2  ^^^2/n  ^2/nj 


At  r  =  a,  =  -p 


With  D  replaced 


r/2 

1/n  / 

f  1 

1 

U/n 

=  -P 


Then 


P 


n 

1/n  / 

/  1 

2 

\^)  B  J 

L2/n 

12 


Substituting  for  C  and  letting  x  =  b/r  and  p  =  b/a  give 


Differentiating  with  respect  to  r  in  equation  (A15)  gives 

WWp2/n_i 


Substitution  into  the  equilibrium  equation  (Al)  yields  the  equation  for 


Og  - 


2/n 


P 


-  1 


(A15) 


(A16) 


From  equation  (A7) 


From  equation  (A9) 


n  p2/n  .  1 


P 


(A17) 


(A18) 
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From  equations  (A4),  (All),  and  (A12) 


(A19) 

2 

^  2 

(A20) 

0  (assumption  (5)) 

(A21) 

T  = 

(A22) 
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^'The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof!’ 

—  National  Aeronautics  and  Space  Act  of  1958 

NASA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 


TECHNICAL  REPORTS:  Scientific  and 
technical  information  considered  important, 
complete,  and  a  lasting  contribution  to  existing 
knowledge. 

TECHNICAL  NOTES:  Information  less  broad 
in  scope  but  nevertheless  of  importance  as  a 
contribution  to  existing  knowledge, 

TECHNICAL  MEMORANDUMS: 
Information  receiving  limited  distribution 
because  of  preliminary  data,  security  classifica¬ 
tion,  or  other  reasons. 

CONTRACTOR  REPORTS:  Scientific  and 
technical  information  generated  under  a  NASA 
contract  or  grant  and  considered  an  important 
contribution  to  existing  knowledge. 


TECHNICAL  TRANSLATIONS:  Information 
published  in  a  foreign  language  considered 
to  merit  NASA  distribution  in  English. 

SPECIAL  PUBLICATIONS:  Information 
derived  from  or  of  value  to  NASA  activities. 
Publications  include  conference  proceedings, 
monographs,  data  compilations,  handbooks, 
sourcebooks,  and  special  bibliographies. 

TECHNOLOGY  UTILIZATION 
PUBLICATIONS:  Information  on  technology 
used  by  NASA  that  may  be  of  particular 
interest  in  commercial  and  other  non-aerospace 
applications.  Publications  include  Tech  Briefs^ 
Technology  Utilization  Reports  and 
Technology  Surveys, 


Details  on  the  availability  of  these  publications  may  be  obtained  from: 

SCIENTIFIC  AND  TECHNICAL  INFORMATION  OFFICE 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

Washington,  D.C.  20546 


